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ABSTRACT 

We propose a method to study the epoch of reionization based on the possible 
observation of 2p-2s fine structure lines from the neutral hydrogen outside the cosmo- 
logical H II regions enveloping QSOs and other ionizing sources in the reionization 
era. We show that for parameters typical of luminous sources observed at z ~ 6.3 
the strength of this signal, which is proportional to the H I fraction, has a brightness 
temperature ~ 20 fiK for a fully neutral medium. The fine structure line from this 
redshift is observable at v ~ 1 GHz and we discuss prospects for the detection with 
several operational and future radio telescopes. We also compute the characteristics 
of this signal from the epoch of recombination: the peak brightness is expected to be 
~ 100 [iK\ this signal appears in the frequency range 5-10 MHz. The signal from the 
recombination era is nearly impossible to detect owing to the extreme brightness of 
the Galactic emission at these frequencies. 

Subject headings: cosmic microwave background — radio lines:general — line: formation — 
radiative transfer 



1. Introduction 

Even though the existence of Hydrogen fine structure lines and their explanation using Dirac's 
atomic theory has been known for close to a century, a Hydrogen fine structure line has never been 
detected from an astrophysical object. An interesting Hydrogen fine structure line is the 2p-2s 
transition. The main difficulty in detecting this line is that the line strength is proportional to the 
population of either the 2p or 2s states which, being excited states, are not so readily populated 
in most astrophysical circumstances. Moreover, the line width of the excited 2p state, which is 
determined by its decay time, is large (99.8 MHz), making the detection of the fine structure line 
a difficult observation. One astrophysical setting where the feasibility of detecting such a line has 
been studied is the interstellar H II regions (see, e.g., Dennison, Turner, & Minter 2005 and refer- 
ences therein; Ershov 1987); in H II regions, the excited levels are populated by recombination. 
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Here we consider two cosmological settings in which the excited levels are populated by either 
recombination or pumping by Lyman-a photons from an external source: (a) The Recombination 
epoch: The Universe makes a transition from a fully ionized to an almost fully neutral medium at 
z — 1089 (Spergel et al. 2006; for details see, e.g., Peebles 1993 and references therein). During 
this era, as the density and temperature of the Universe drops, recombination is stalled owing to a 
high Lyman-a radiation density and progresses either by the depopulation of the 2p state owing to 
redshifting of the photons out of the line width or the 2-photon decay of the 2s state. This results 
in a significant 2 s and 2p level population during the recombination era. (b) The Reionization 
epoch: Recent observations suggest that the universe made a transition from nearly fully neutral to 
fully ionized within the redshift range 6 £ z ^ 15 (Page et al. 2006; White et al. 2003; Fan et al. 
2002; Djorgovski et al. 2001; Becker et al. 2001). It is widely believed that this 'reionization' 
was achieved by the percolation of individual H II regions around the sources of reionization. The 
nature of these sources is not well understood: they might be Pop ill stars, active galactic nuclei 
or star-forming galaxies. During this epoch, a signal from the 2p — 2s fine- structure transition 
might originate from either within the cosmological H II regions or from the almost fully neutral 
medium surrounding the H II region. The level population of the first excited state in the former 
case would be largely determined by recombinations and in the latter case by Lyman-a photons 
from the central source. We shall show below that for most cases of interest the fine-structure line 
from within the cosmological H II region might be negligible as compared to the signal from the 
regions immediately surrounding the H II region. 

Throughout this work we adopt the currently-favoured ACDM model: spatially flat with 
Vt m = 0.3 and tt A = 0.7 (Spergel et al. 2006; Riess et al. 2004; Perlmutter et al. 1999) with 
Q b h 2 = 0.022 (Spergel et al. 2006; Tytler et al. 2000) and h = 0.7 (Freedman et al. 2001). 



2. Fine-structure lines from the reionization epoch 

Subsequent to the recombination of the primeval baryon gas at redshift z ~ 1089 (Spergel 
et al. 2006) and the transformation of the gas to an almost completely neutral state, it is believed 
that the gas was reionized during epochs corresponding to the redshift range 6 £ z £ 15. WMAP 
measurements of cosmic microwave background radiation (CMBR) anisotropy in total intensity 
and polarization have been used to infer that the baryons were likely neutral at redshifts z £ 12 — 
15; however, the detection of CMB polarization anisotropy requires substantial ionization by about 
z ~ 11 (Page et al. 2006). Observationally, the Gunn-Peterson (GP) test shows that the universe is 
highly ionized at redshifts lower than z ~ 5.5; the detection of GP absorption at greater redshifts 
suggests that the neutral fraction of the intergalactic hydrogen gas rises to at least 10~ 3 in the 
redshift range 5.5 £ z £ 6, and that reionization was not complete till about z ~ 6 (White et al. 
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2003; Fan et al. 2002; Djorgovski et al. 2001; Becker et al. 2001). However, from the GP test 
alone it is not possible to infer the neutral fraction of the medium; it only gives a rather weak lower 
limit of ~ 1CT 3 on the neutral fraction of the universe for z £ 6. From other considerations it is 
possible to put more stringent bounds on the neutral fraction; for example, Wyithe & Loeb (2004) 
obtain a lower limit of 0.1 on the neutral fraction of the universe at z ~ 6.3 (see also Mesinger & 
Haiman 2004). 

Our understanding of the nature of the sources that caused the reionization is far from com- 
plete. The transition from an almost completely neutral gas to a highly ionized gas during redshifts 
is a key problem in modern cosmology and considerable theoretical and experimental 
efforts are currently directed at this unsolved problem. Here we propose a new method, based 
on the 2p — 2 s fine- structure transition, to determine the evolution of the neutral fraction of the 
intergalactic medium within this epoch. 

Owing to fine structure splitting, the two possible transitions between the 2s and the 2p states 
are: 2p 1 / 2 -2s 1 / 2 at a frequency ~ 1058 MHz that has an Einstein A coefficient 1.6 x 10~ 9 s^ 1 and 
2p 3 / 2 -2si/ 2 at a frequency u ps ~ 9911MHz that has an Einstein A coefficient 8.78 x 10~ 7 s _1 . 
The Einstein A coefficient for the latter transition is more than an order of magnitude greater than 
the former; therefore, in this work we consider only the 2p 3 / 2 -2si/ 2 transition and hereinafter we 
refer to this specific transition simply as the 2p-2s transition. 

The ionizing UV photons from sources in the reionization era create 'Stromgren spheres'. 
Whereas the gas in the cosmological H II regions are highly ionized by the photons, the ioniza- 
tion level of the gas beyond the Stromgren spheres is determined by the history of the gas, the 
density, and the mean specific intensity of the background ionizing photons, which includes both 
the photons diffusing out of the Stromgren spheres as well as the background radiation field. In 
this work we assume that the ionizing sources at these high redshifts are AGNs and in illustrative 
examples adopt parameters of a few QSOs that have been observed at z ~ 6. The photons at 
the Lyman-a transition frequency (here and throughout, unless otherwise specified, we shall con- 
tinue to refer to frequencies between Lyman-a and Lyman-limit as 'Lyman-a) from a high redshift 
QSO escape the mostly-ionized Stromgren sphere and are strongly scattered and absorbed in the 
medium beyond. The population of the 2p level in this region is determined by (a) the intensity of 
Lyman-a photons from the central source, (b) recombination rate of free electrons, (c) absorption 
of CMBR photons by electrons in the 2s state (it is assumed here and throughout this work that the 
only radio source at high redshifts is the CMBR) and (d) collisional transition from the 2s state. 
The 2s state is populated via (a) the recombination rate of free electrons, (b) collisional transfer 
of atoms from the 2p state, (c) the spontaneous decay of the 2p state, and (d) transition from the 
2p state stimulated by CMBR photons. Additionally, the absorption of photons from the central 
sources, with energy equal to or in excess of the Lyman-/? transition, would result in electronic 
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transitions to the second excited state, which could be followed by spontaneous decay to the 2s 
state. (It might be pointed out here that both the 2s and 2p states could also be populated by atoms 
cascading from excited states with n > 3. In particular, all photons absorbed from Is states to any 
excited state can directly de-excite to the 2s level. However, the rate of transition from Is to any 
excited state is roughly oc l/n 3 (e.g. Rybicki & Lightman 1979) and, therefore, we include only 
the most dominant transition in each case.) 

The population of the ground state is denoted by n ls . We denote the level populations of the 
the two states — 2s\i 2 and 2p 3 / 2 — by the number densities n 2s and n 2p ; these may be solved for, 
respectively, from the following two equations of detailed balance: 

fa B rii + cB 2p2s n 2p n CMBIi + C ps n,jn 2p + A 2p2s n 2p + cn ls p 32 J B 13t p<j> 13 (v)n a (v)di/ 

= A 2sls n 2s + C sp riin 2s + cB 2s2p n 2 

S^CMBR) (1) 

(1 - 2/)a B n 2 + cB 2s2p n 2s n CMBR + C sp riin 2s + cn ls J B l2 ^ a (j) l2 (v)n a (ij)dv 

= A 2p \ s n 2p + C ps nifi 2p + cB 2p2s n 2p n C MBK- (2) 

Here / is the fraction of all the atoms that recombine to the 2s state. In equilibrium f — 1/3 
as the n = 2 state splits into three doublets: 2pi/ 2 , 2p 3 i 2 and 2s\i 2 . olb is the recombination 
coefficient and n» is the density of the ionized gas. B 2p2s = B 2s2p = c 2 / (8irv ps ) A 2p2s is the Einstein 
B coefficient for the 2p 3 / 2 -2si/ 2 transition in terms of the corresponding Einstein A coefficient 
A 2p2s (note that the two B coefficients are equal as the two states have the same degeneracy). 
Cps — C S p — 5.31 x 10 4 cm 3 s 1 is the rate coefficient of transition owing to collisions with 
electrons. n CMBH is the number density of CMBR photons within the transition line width. n a (v) 
is the number density (per unit frequency) of photons with frequency equal to or larger than the 
Lyman-a frequency (and smaller than the Lyman-limit frequency) at any location; 13 and <j) 12 
are, respectively, the line profiles corresponding to the Lyman-/? and Lyman-o; transitions, and 
p 32 is the probability for the electron transition to the 2 s state following excitation to n = 3 via 
absorption of a Lyman-/? photon. We have not included the induced Lyman-a transition because 
the number density of atoms in the 2p state is negligible as compared to that in the Is state. p 32 is 
the probability that an atom in the third excited state (3p) will decay to the 2s state. A 2sls is the 
Einstein A coefficient corresponding to the 2-photon decay of the 2s state. Other symbols have 
their usual meanings. 

Owing to the fact that the mostly neutral medium in the vicinity of the cosmological Strom- 
gren spheres is optically thick to Lyman-a scattering, the Lyman-a photons from the decay of 
the 2p state are strongly scattered by the gas. Therefore, n a (u) will contain contributions from 
both the Lyman-a photons from the central source as well as the Lyman-a photons that arise 
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from recombinations outside the Stromgren sphere and are multiply-scattered therein: n a (u) = 
n source^^ _|_ ^rec^^ ne gi ec t the multiply- scattered Lyman-a photons from the Stromgren 
sphere that have been reprocessed via recombination within the Stromgren sphere because these 
would be redshifted redward of the Lyman-a line before encountering the boundary of the Strom- 
gren sphere. The scattering of recombination photons in an optically thick, expanding medium is a 
complex problem (Field 1959; Rybicki & Dell' Antonio 1994). One of its first applications was to 
study the recombination of primeval plasma (Peebles 1968; Zeldovich, Kurt & Sunyaev 1969). In 
these analyzes it was implicitly assumed that apart from 2-photon decay, in an expanding universe 
the dominant effect that results in resonant photons ceasing interaction with the gas, and leaving 
the system, is its redshifting out of the line profile. The effect of scattering off the moving atoms 
was deemed to be either negligible or at best comparable. This assumption has been borne out by 
more recent detailed analysis that have taken into account the effect of scattering on the photon 
escape (Krolik 1990). Taking only the redshift as the main agent of photon escape, it can be shown 
that the net effect of the scattering of a resonance photon before it drops out of consideration is to 
reduce the decay time of the 2p state from A 2 \ = 6.2 x 10 8 s _1 to A 2 i/t gp (Zeldovich et al. 1969; 
for more recent work see, e.g., Chluba, Rubino-Martin & Sunyaev 2007, Seager, Sasselov & Scott 
1999 and references therein; we give a concise derivation in Appendix A). In this expression, the 
Gunn-Peterson optical depth r GP = [3/(8irH)]A 2p i s \ 3 a ni s . 

Similar complications exist in computing p 32 , the probability that an atom in the third excited 
state will decay spontaneously to the 2s state, in an optically thick medium. In optically thin media, 
P32 = A 32 / (A 32 + A 3 i); on the other hand, in an optically thick medium, we would need to take 
into account the 'trapping' of the Lyman-/? photon owing to resonant scattering. The effect of this 
scattering in an optically thick medium would be that the fraction of photons that decay directly 
to the ground state are reabsorbed 'locally' to the third excited state and, therefore, all photons 
absorbed to the third excited state result in an Ha photon and an atom in the 2s state. This means 
that the appropriate value of p 32 is close to unity in an optically thick medium: in this work we 
assume p 32 = 1. 

The astrophysical setting in which we seek solutions to the algebraic equations above is cos- 
mological H II regions at high redshift. In particular, we are interested in the signal from the neutral 
region surrounding the cosmological H II region. For a fully neutral inter-galactic medium (IGM) 
at z ~ 6.5, r GP ~ 6 x 10 5 . If we adopt spectral luminosities corresponding to QSOs observed at 
these high redshifts, it may be shown that the populating of the 2p state via direct recombinations 
from the free-free state, collisional transfer from the 2 s state, and upward transitions from the 2s 
to 2p state arising from absorption of background CMBR photons may all be neglected. (The rele- 
vant parameters are: n,i ~ n b ~ 2.8(1 + z) 3 cm~ 3 in the H II region surrounding the sources, with 
rii expected to be much smaller in the neighboring mostly neutral medium; the number density 
of CMBR photons that might cause a 2p-2s transition is ~ 5(1 + z) cm -3 ; the number density 
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of Lyman-a photons from the central source, assuming luminosities typical of SDSS quasars at 
z ~ 6.5 (more details in §4), is n a ~ 10~ 4 cm -3 . First, for these parameters, the dominant process 
that populates the excited state is the pumping by Lyman-a photons. Second, it may be readily 
verified that for these plausible values for the parameters the signal expected from the H II region 
surrounding the central source is negligible as compared to the signal from the surrounding neu- 
tral region.) Given that the Lyman-a flux from the central QSO is the dominant causative factor 
for populating the 2p state, the two equations CD and [2]) that determine the level populations is 
simplified. The number density of atoms in the 2p state is given approximately by: 

n 2p - n ls T a T GP /A 2 p ls , (3) 

where n ls = f ncu n h , with / neu denoting the neutral fraction and n b ~ 2.7 x 1CT 7 (1 + z) 3 cm~ 3 
is the number density of baryons in the IGM. T a = J B 2p i s cf)i 2 {y)n s ° nrce {v) (in units of s" 1 ) is 
the transition rate to the 2p state owing to the Lyman-a photons from the central source, where 
n s ° urce (v) is the number density of Lyman-a photons from the central source alone. 

Similarly, the dominant process that determines the population of the 2s state is the absorption 
of Lyman-/? photons (for details, see the discussion above and §4) and the subsequent decay to the 
2s state. The number density of atoms in the 2s state is: 

n 2s ^ n ls T f3 /A 2sls , (4) 

where Tp = f i?3pi s 0i3(^)n^ ource (z/) (in units of s -1 ) is the transition rate to the 2s state owing to 
the Lyman-/? photons from the central source. 



3. Fine structure lines from the epoch of recombination 

The universe made a transition from fully ionized to nearly fully neutral at z ~ 1089 (Spergel 
et al. 2006; Peebles 1968; Zeldovich et al. 1969). This transition is mainly accomplished by the 
2-photon decay of the 2s state and the slow redshifting of the Lyman-a photons which deplete the 
2p state (see Peebles 1993 and references therein for a detailed discussion). The Saha ionization 
formula, valid for thermodynamic equilibrium conditions between the hydrogen level populations 
and free electrons, is a poor approximation for studying the epoch of recombination. A good ap- 
proximation for studying this transition is to assume that all states, excepting Is, are in equilibrium 
with the CMBR (matter temperature to a very good approximation remains equal to the CMBR 
temperature throughout this transition) (Seager et al. 1999; Peebles 1968). However, in this ap- 
proximation, where the matter temperature and all transitions, excepting the Lyman-a line, are in 
thermal equilibrium, the 2p-2s signal is unobservable because the excitation temperature for this 
transition equals the background radiation temperature. Even though this second approximation 
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might be useful for studying the evolution in ionization fraction, it will be strictly true only if the 
dominant mechanisms that determine the level populations of the 2p and the 2s states are either 
interaction with the CMBR photons or collisions between atoms. As there are a variety of other 
processes relevant to the determination of the level populations — for example, the free decay of 
either of the two states — a small deviation from equilibrium is expected in the 2s and 2p level 
populations and it is our aim here to compute it. 

One approach to this problem is to simultaneously solve for the level populations of the 2p 
and the 2s states as well as the change in the ionization fraction. However, assuming thermal 
equilibrium between these two states is a good approximation for solving the evolution of ioniza- 
tion. Therefore, the approach we have adopted is to solve for the evolution of ionization using the 
method of Peebles (1968), and use the resulting ionized/neutral fraction to solve for the populations 
of the 2s and the 2p states using detailed balance. The resulting equations are: 

fa B n? + cB 2p 2sn 2p n CMBR + C ps nin 2p + A 2 p2 S n 2p + n ls A 2sU exp(-(B 1 - B 2 )/(k B T C MBR)) 

= A 2s i s n2 S + C sp nin 2s + cB2 S 2 P n2 S n CMBR + [3 c n 2s (5) 

and 

(I - 2f)a B rf + cB2 S 2 P n2 S n CMBR + C sp nin 2s 

— A 2 pl s n2p/ T GP + C vs n(ll2 v + cB2p2sTl2 p TlcMBR + Pc n 2p(6) 

where B x = 13.6 eV and B 2 = 2.4 eV are, respectively, the ionization potentials of the ground 
and the first excited states, k B is the Boltzmann's constant, and T C mbr is the temperature of the 
CMBR. p c is the rate at which the CMBR photons cause a bound-free transition of electrons from 
the n = 2 state (2s or 2p states). The various other terms in these equations have the same meanings 
as in the previous case. The main difference is that the CMBR photons and baryonic matter at the 
time of recombination are hot and dense enough to directly affect the level populations of the 
excited states by ionizing the excited state, and the two-photon capture to the excited state is not 
completely negligible (see Peebles 1993 and 1968 for details of the different physical processes 
that are relevant at this epoch). 



4. Expected brightness of the fine structure lines 

The brightness temperature in the 2p-2s transition is: 

\rp rp rp n n 2p (0) h p A 2p 2s ^ e ( 1 + zf C 

IXlb = J-b — -L CMBR — 9* —, — 777 s (,1 — JCMBR/ J- ex) = ^ex\J- ex ~ ^CMBrJ {') 

01TK B li [Z ) 

Here n 2p (0) = n 2p /(l + z) 3 and A e is the rest wavelength of the 2p-2s transition. T ex = 
(h p u e / k B )[n 2p / (n 2s — n 2p )] is the excitation temperature corresponding to the transition, where 
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h p is the Planck constant. r ex is the optical depth of the source in the 2p-2s signal, takes into 
account the selection rules for transitions between the 2p 3 / 2 and 2si/ 2 levels. Given the selec- 
tion rules there are three allowed transitions between these two states; they occur at frequencies 
v ~ 9852, 9875, and 10029 MHz (see, e.g., Ershov 1987). The first two transitions are blended by 
the natural width of the line, which is approximately 100 MHz, but the third should be observable 
as a distinct line. This implies that g+ = 2/3 if the observing frequency is ~ 9900/ (1 + z) MHz. 



4.1. Expectations for signals from cosmological H n regions 

For computing the strength of this signal in a typical case, we adopt observed parameters 
of QSO SDSS J1030+0524 (see, e.g., Wyithe & Loeb 2004 and references therein), which is at 
redshift z = 6.28 and shows no detectable flux beyond the QSO Stromgren sphere: the radius 
of the Stromgren sphere has been estimated to be R ~ 4.5 Mpc (Mesinger & Haiman (2004) 
argue that the size of the Stromgren sphere could be roughly 30% higher; this makes no essential 
difference to our results). It may be noted here that the QSO SDSS Jl 148+5251 also has similar 
parameters. We arrive at an estimate of n^ ource (z/) by assuming that L a photons per second, with 
wavelengths corresponding to the Lyman-a transition, are emitted by the central source in an 
effective frequency range Az/ sourcc , and that the photons are absorbed at a radial distance R. This 
leads to: 

< urcc = L a /(A7rR 2 c) x - J—. (8) 

^ ^source 

We assume typical values: L a = 10 58 s _1 , R = 4.5 Mpc and Az/ source = 30, 000 km s -1 . These 
lead to the following estimate for the transition rate, r Q , in the gas at the boundary of the Stromgren 
sphere arising due to the photons from the central QSO: 

r a ~8x 10~ 10 s _1 . (9) 

It should be noted that the mean specific intensity of Lyman-a photons in the IGM would also 
give a non-zero signal. Assuming a mean specific intensity of ~ 10~ 21 erg cm -2 sec -1 Hz -1 sr _1 
(this might be needed to couple the HI spin temperature to the matter temperature; see e.g. Madau, 
Meiksin, & Rees 1997), the expected signal is many orders of magnitude smaller than we have 
computed from the outskirt of bright sources. If we assume that the central sources are continuum 
emitters in the frequency range between Lyman-a and Lyman-limit frequencies, we may assume 
similar parameter values for computing the expectations for Tp. In that case, 
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where fp and f a are the oscillator strengths of the Lyman-/3 and Lyman-a transitions, respectively. 
From Eqs. ©, © and (fT0|): 

n 2p f a A 2s 

For a completely neutral medium at z ~ 6.4, r GP ~ 6 x 10 5 , which may be used to show that 
n 2s n 2 P f° r the parameters of IGM in the redshift range of interest: 6.4 <, z ^ 10. This implies 
that in the outskirts of the Stromgren sphere surrounding the QSO, the transition is expected to 
be observable as an absorption feature against the CMBR. Using Eqs. © and © in Eq. ©, the 
observable brightness temperature is estimated to be: 

AT B ~ -20 iiK ; (12) 

where / neu is the neutral fraction of hydrogen outside the Stromgren sphere, and might be close to 
unity. 

We have adopted parameters typical of QSOs observed at the edge of the reionization epoch 
in estimating the above temperature decrement. The main uncertainty above is in the estimation of 
the 'Lyman-a' flux from the central source, and as the observed temperature decrement is directly 
proportional to this flux from the central source, this constitutes a major uncertainty in reliably 
computing the expectations for the signal. For QSOs that have strong Lyman-a and Lyman-/? 
lines, the blueward side of the lines will be strongly absorbed in the medium just beyond the 
Stromgren sphere (and this has been observed to happen in many cases), provided that the blueward 
side photons have not been redshifted to frequencies smaller than the Lyman-a frequency while 
transversing the Stromgren sphere. In the case of QSOs that have large line fluxes and small 
Stromgren spheres, the Lyman-ct luminosity L a might be underestimated. 



4.2. Expectations for the fine-structure line from recombination 

Using equations © and © the level populations of the 2p and the 2s states may be computed. 
Solving for the level populations and using equation ©, we have computed the expected signal 
from the recombination epoch; the expected signal is shown in Fig. 1. The fine structure line 
transition is expected to be an absorption feature, with a maximum temperature decrement of order 
100 /iK at an observing frequency of about 7 MHz. The width of the decrement in frequency space 
corresponds to a redshift span Az ~ 200, which is roughly the width of the visibility function at 
recombination. 
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5. Prospects for the detection of the cosmological fine structure lines 

As Eq. (fT2l) shows, a detection of the fine structure line in the outer regions of cosmological 
H II regions is potentially a probe of the cosmological neutral hydrogen density in the vicinity of 
QSOs in the reionization epoch, and might be a tool for the investigation of the evolution in the 
neutral fraction with cosmic epoch through the reionization era. Given the cosmological impor- 
tance of such a measurement, and the fact that there does not exist many reliable methods for the 
detection of H I at high redshifts (see, e.g., Barkana & Loeb 2001), the detection of the 2p-2s line 
transition in the cosmological context assumes additional importance. 

In deriving equation © the Hubble expansion was assumed to be the only cause for the 
velocity width in the observed line. However, an important contribution to the velocity dispersion 
in the line in this case is the natural width of the fine structure line: owing to the rapid decay of 
the 2p state, the natural Lorentzian width in the rest frame of the gas is Vior — 100 MHz (see, 
e.g., Dennison et al. 2005). Therefore, the peak brightness in the observed line profile might 
be suppressed by a factor V cxp /Vi, OI , where V exp is the line-of-sight peculiar velocity dispersion 
owing to the Hubble flow across the region being observed. However, in the case V cxp £ VLor this 
suppression has a negligible effect. For QSOs at redshift z ~ 6.5, the natural Lorentzian width 
of the fine structure line is equivalent to the peculiar Hubble flow across a proper line-of-sight 
distance of ~ 4 Mpc. This distance is approximately the size of the Stromgren spheres around 
QSOs at that redshift; therefore, the natural width of the line does not significantly diminish the 
expectations, given by equation ©, for the peak brightness temperature. A second inference is 
that in the case of a QSO that has a smaller Stromgren sphere, the increase in mean brightness 
temperature is roughly proportion to the inverse of the radius of the Stromgren sphere: 1/R, and 
not 1/R 2 . 

Another assumption that was made while deriving Eq. © is that the only radio frequency 
radiation that needs to be considered for the determination of the level populations and brightness 
temperature decrement is the CMBR. It may be that the central ionizing source, which may be a 
QSO, is radio loud. There may also exist radio sources behind the observed Stromgren sphere and 
within the angular region over which the Lyman-a flux from the QSO is appreciable. Equation 
CD) may be modified, to account for this, by replacing T C mbr with Tcmbr + T B , where T B is the 
brightness temperature of the radio source at wavelengths corresponding to the rest frequency of 
the transition, which is ~ 9 GHz. The result would be to enhance the brightness temperature of 
the line signal by ~ T b /T C mbr- If the radio source is unresolved, it is appropriate to estimate the 
expected signal in terms of the optical depth. The optical depth corresponding to equation (fl~2l) is 
~ 10~ 6 . 

We now discuss the feasibility of the detection of the fine structure line towards SDSS J1030+0524. 
We shall assume that the neutral fraction, / ncu , outside the Stromgren sphere of this QSO is unity, 
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consistent with the measured GP trough. The redshifted fine structure line would be expected 
at ~ 1.36 GHz. The observed line 'width', considering natural broadening and the Hubble flow 
across the Stromgren sphere, is expected to be approximately Ais ~ 100MHz/(l + z); using 
z = 6.28, we obtain Ais = 13.7 MHz. The fine structure line would be expected to originate in a 
shell that is roughly the size of the Stromgren sphere, and fall of as 1/r 2 beyond this shell, where r 
is the distance from the QSO. In the case of SDSS J1030+0524, the angular size of the Stromgren 
sphere is expected to be 15 arcmin. 

The frequency of the redshifted line is in the observing bands, and the line width is within the 
spectral line capabilities, of several currently operational telescopes. However, large-collecting- 
area arrays like the Giant Metrewave Radio Telescope (GMRT) have large aperture antennas of 
45-m diameter and, therefore, poor surface brightness sensitivity for such extended structures. 
The Australia Telescope Compact Array (ATCA), with 22-m antennas, has reasonable brightness 
sensitivity for this problem and, additionally, operates with 128-MHz bandwidths in spectral line 
mode. At v ~ 1.3 GHz, the ATCA has a system temperature T sys ~ 25 K and antenna efficient 
K — 0.1 K Jy -1 . The five movable 22-m diameter antennas may be configured into an ultra- 
compact 2-D close-packed H75 (75-m maximum baseline) array, and this would yield a number 
of baselines sensitive to the 15-arcmin scale fine-structure line signal. The brightness sensitivity 
of this array, for a 8-arcmin scale structure, is 400 /iK in 6 hr integration time. The brightness 
sensitivity in Fourier synthesis images could be enhanced somewhat, by factors of a few, by ap- 
propriately weighting the baselines to match the synthetic beam to the expected structure scale; 
however, the required integration times for a detection are still in the ball park of 10 3 hr. 

The detection of the signal from cosmological H II regions, however, might be feasible using 
facilities under construction or planned for the near future, like the xNTD in Australia or the Square 
Kilometer Array (SKA). These arrays would have smaller antenna sizes — making the detection of 
these large-angular-scale structures detectable in interferometers — and significantly more numbers 
of antennas, giving more numbers of short baselines that would usefully respond to the large- 
angular-scale fine-structure line. However, the array configuration designs would have to factor in 
the extraordinary brightness temperature sensitivity requirements for this demanding observation. 

As an example, the SKA might consist of about 10 4 12-15 m class reflector antennas, with 
aperture efficiency of 60%, and the system temperature at 1.4 GHz might be about 20 K. Assuming 
that the visibilities are optimally weighted so that the synthetic beam of the Fourier synthesis array 
is matched to the source size of about 15 arcmin FWHM, the line strength would have a peak of 
about 20 //Jy. A 5-er detection of the fine structure line towards a typical QSO, in a reasonable 
integration time of about 1 hr, will require that about 250 baselines (just 0.001% of the total) be 
within about 50 m. 

The signal from the recombination epoch is observable in the frequency band 6-8 MHz as 
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a broad decrement in the brightness temperature of the extragalactic background sky, and would 
be extremely small compared to the orders of magnitude more intense Galactic non-thermal emis- 
sion as well as the average low-frequency background brightness temperature arising from the 
numerous extragalactic radio sources. This decrement may be considered to be a distortion to the 
CMBR spectrum at long wavelengths, and would be an all-sky cosmological signal. However, 
owing to the ionosphere, the frequency range in which this feature is expected to appear is too 
low a radio frequency to be easily accessible using ground-based observatories. Therefore, even 
though the observation of this signal might be yet another tool to probe the epoch of recombi- 
nation, new custom-made instruments, which should presumably operate from space and above 
the ionosphere, will require to be built if a detection of this signal is to be attempted. An addi- 
tional cause for concern is that the Galactic and Extragalactic background radiations might have 
low-frequency spectral turnovers at these frequencies, as a result of free-free absorption as well as 
synchrotron self-absorption, and these would result in significant spectral features in the band that 
would require a careful modelling in order to detect any CMBR decrement feature arising from 
fine structure transition absorption. Interference from terrestrial man-made transmitters, as well as 
from auroral phenomena and solar system objects would also be an issue. 



6. Summary and Discussion 

We have discussed the possibility of detecting the hitherto undetected fine structure line of 
2p-2s transition in two cosmological settings: the epoch of recombination at z ~ 1089 and the 
epoch of reionization at z ~6-15. 

The expectations for the line from the environments of ionization sources in the epoch of 
reionization are interesting and worthy of attention as a novel tool for the investigation of the 
reionization process and the cosmological evolution of the gas. The signal is expected to be ob- 
servable as an extended and weakly absorbing source, which causes a decrement in the brightness 
of the background CMBR, and may be detected by interferometers as a negative source akin to 
the Sunyaev-Zeldovich decrements observed along the lines of sight through hot gas in clusters 
of galaxies. Detection of the 2p-2s line signal from the outskirts of cosmological H II regions at 
different redshifts within the reionization era may serve to determine the neutral fraction of the 
medium during the epoch of reionization, which is a quantity of significant interest in modern cos- 
mology. In particular, we have computed a representative signal strength by adopting parameters 
typical of a QSO observed at z ~ 6. These QSOs show GP troughs in their spectra; however, 
owing to the weakness of the GP test, the spectra have only been useful in setting weak limits 
on the neutral fraction (constraining the neutral fraction to be £ 1CT 3 ) outside the observed H II 
region. We have shown that for a fully neutral medium the line peak may reach ~ 20 jiK, which is 
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potentially observable by radio interferometers that are being designed today. 

Other interesting probes of the reionization epoch include detecting 01 line from this epoch 
(Hernandez-Monteagudo et al. 2006). 

It is of interest to compare the relative difficulty associated with detecting the neutral hy- 
drogen in this indirect way, using the redshifted fine structure line, with direct imaging of the 
redshifted 21 -cm line from neutral hydrogen during the epoch of reionization (see, e.g., Sethi 
2005; Zaldarriaga, Furlanetto, & Hernquist 2004). The 'all-sky' H I signal might be detectable 
with a peak strength of ~ 50 mK; however, it could be very difficult to detect owing to calibra- 
tion and foreground contamination issues (e.g., Zaldarriaga et al. 2004; Shaver et al. 1999, and 
references therein). A better approach might be to attempt to detect the fluctuating component 
of the sky signal, which could have peak intensities of ~ a few mK at observing frequencies of 
v ~ 100-200 MHz (Zaldarriga et al. 2004; Shaver et al. 1999). This translates to roughly the 
same signal strength (specific intensity) as we have obtained in Eq. (fT2l) for the fine structure line. 
This is not entirely unexpected: even though the level populations of the excited states are much 
smaller as compared to the ground level population needed in computing the HI signal, the A coef- 
ficient of the fine structure transition we consider here is roughly 8 orders of magnitude larger than 
the HI hyperfine transition A coefficient. Currently, there are many ongoing and planned radio 
interferometer experiments for the detection of the redshifted H I emission/absorption from the 
epoch of reionization {e.g., Pen, Wu, & Peterson 2004; the LOFAR project at www . lof ar . org; 
the MWA project at www .hay stack, mi t.edu/ast/ arrays /mwa/ site/ index . html). 
The detection of the H I signal is firstly more difficult — requiring greater sensitivity — because the 
typical frequency width of the signal ~ 0.5 MHz (Zaldarriga et al. 2004), which is far smaller than 
the typical width expected for the fine structure line (~ 10 MHz). Second, the fine structure line 
is expected at 1 .4 GHz, where the sky background temperatures are significantly lower than in the 
100-200 MHz band, making the telescope system temperatures lower. An additional advantage 
of the the indirect detection is that the redshifted fine structure line appears at higher frequencies 
( ^ 1 GHz), which are relativity free of interference as compared to the low frequency band of 
100-200 MHz. The main advantage of the direct detection is that unlike the method we suggest 
here, it is independent of the existence of strong Lyman-a emitters at high redshifts. Another ad- 
vantage of direct detection is that it may be detected 'statistically' and such a detection might be 
achieved with greater ease than direct 'imaging' (e.g Zaldarriaga et al. 2004; Bharadwaj & Sethi 
2001); however, the foreground subtraction problem becomes a severe constraint for a statistical 
detection. To summarize: if strong Lyman-a emitting sources are present at high redshifts, they 
would facilitate the indirect detection of the neutral hydrogen via enabling the detectability of the 
fine structure line. The imaging issues and problems associated with the detection of this signal 
appears to be less of a challenge as comparison to the direct detection of redshifted H I from those 
epochs. 
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Another astrophysical context in which the fine structure line might be detectable is the envi- 
ronments of high redshift galaxies, which are strong Lyman-a emitters. As equation (PT2l) shows, 
the observed signal depends on the neutral fraction outside the Stromgren sphere. Therefore, de- 
tection of this signal would constitute an alternate probe of the neutral fraction of the IGM at large 
redshifts. 

The aim of this work has been to examine the detectability of the fine structure line in cos- 
mological contexts, to point out the cosmological significance of detections, and spawn work that 
may refine the modelling presented herein and improve the case for appropriate design of future 
telescopes, which might enable the detection of the fine structure line towards multiple sources in 
the reionization era. 



Appendix A: Photon distribution function 

The evolution of the photon distribution function, neglecting the effect of scattering off mov- 
ing atoms, is (see e.g. Rybicki & Dell' Antonio 1994): 

- vH-^j = A 2pls n 2 p(j) u - cB u n ls (j) v n u . (13) 

Here B v = 3/(87r)c 2 /^A 2p i s and H = a/a. The equation may be written as: 

^ _ vH ^L = -1/ T (n u - nA (14) 
dt dv 1 y v ' y ' 

Here r = 1/ (cB u (p u ni s ) and /* = A 2p i s ri2p/ (B u ni s c). The equation above lends itself to a ready 
interpretation. If the second term on the left hand side (which is owing to the expansion of the 
Universe) was absent, the distribution function will approach /* on a time scale ~ r, where r ~ 
3 x m s s <C a/a ( n ls here has units cm -3 ). It may be readily verified that for the recombination 
epoch and also the epoch of reionization, r <C l/H, excepting when the neutral fraction of the 
medium is very small. We work here with the assumption that the neutral fraction is always large 
enough so that r <C l/H. In this case the solution to Eq. (fT4l) may be simplified: to leading order 
the distribution function approaches /* and the first order term (of the order of tH) represents the 
slow time variation of the distribution function owing to the expansion of the Universe. In this 
approximation, one may write the solution to Eq. (fT4l) as: 

n v ~ n* - rH(t)n*. (15) 

Using these equations we may proceed to prove the contention that the net effect of the 'trapping' 
of Lyman-a photons is to reduce the decay time of the 2p state by a factor t gp . The ls-2p transition 
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rate, which is given by Eq. ©, may be solved using Eq. (fT5l) : 

cn ls J B v §{y)n<y)dv ~ A 2pls n 2p - A 2pls n 2p /r GP . (16) 

The first term on the right hand side cancels with the decay term of the 2p state on the right hand 
side of Eq. ©, and, therefore, the net effect of the scattering of recombination photons is to reduce 
the decay time of the 2p state by a factor r GP . It may be pointed out that the condition needed to 
derive the above expression roughly translates to the condition that r GP 3> 1. For the reionization 
case, this requires that the neutral fraction ^ 10 5 . In the case of primordial recombination, it 
leads to an even weaker condition that the neutral fraction is ,> 1CT 7 . 
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Fig. 1. — The expected brightness temperature decrement in the background radiation, owing 
to the fine structure transition in gas at the recombination epoch, is plotted versus the observing 
frequency. 



